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In heavy quarkonium production, the measured ratio Rcc = o-[J/tp + cc + X]/a[J/tl} + X] at 
B factories is much larger than existing theoretical predictions. To clarify this discrepancy, in 
nonrelativistic QCD (NRQCD) we find the next-to-leading-order (NLO) QCD correction to e'^e~ — > 
J/^ + 99 can enhance the cross section by about 20%. Together with the calculated NLO result 
for e^e~ J /ijj -\- cc, we show that the NLO corrections can significantly improve the fit to the 
ratio Rcc- The effects of leading logarithm resummation near the end point on the J/i/> momentum 
distribution and total cross section are also considered. Comparison of the calculated cross section 
for e^e~ — > J/^j) + gg with observed cross section for e"'"e~ J/tp + non{cc) is expected to provide 
unique information on the issue of color-octet contributions. 

PACS numbers: 13.66.Bc, 12.38.Bx, 14.40.Gx 



In recent years a number of challenging problems in 
heavy quarkonium production have appeared[lj. Aside 
from the J/^ production cross sections and polariza- 
tions in hadron-coUisions at the Tevatron, charmonium 
production in e+e^ annihilation at B factories^, [sj 
also conflicted with theoretical predictions. The ob- 
served double charmonium production cross sections for 
e'^e'' — !■ J / i^ric{xco) were larger than the LO calculations 
in NRQCD [4] by an order of magnitude [5], and later it 
was found that these discrepancies could be largely re- 
solved by the NLO QCD corrections (see 0] for J/it>ric 
and for J/ipXco) with relativistic corrections 0, 
For the J/ip production associated with an open charm 
pair e+e^ J/ip + cc, the NLO QCD correction [ll| was 
also found to significantly enhance the cross section (see 
also [l^). and reduce the large gap between experiment 
and the LO calculations 13 1. 

Another important issue concerns the ratio 



Rcc — 



■X] 



J/^ + X] 



(1) 



Belle found first 0.59tEJ i^±O.120, and later R^^ 

0.82 ± 0.15 ± 0.14 [14^. On the contrary LO NRQCD jll, 
15 | and light-cone pQCD predictions[ld| for the ratio are 
only about 0.1 — 0.3. The color evaporation model gives 
a value of only 0.06[l7|. 

In NRQCD, (t[J/V' + X] includes color-singlet contri- 



butions (T[J/i/;(3S'pl) 



and (t[J/i/;(3S'W) -t- gg], and 



color-octet contributions a[J/4'{'^P 



J ' 



'S^^^) + g]. 



Con- 



tributions of other Fock states are suppressed by either 
as, the strong coupling constant, or v, the relative veloc- 
ity between quark and antiquark in heavy quarkonium. 



was calculated at LO in a.hi 



and an apparent enhancement at the J /ip maximum en- 
ergy was predicted. But experiments did not show any 
enhancement at the end point. The resummations of the 



V expansion and log(l — z) where z — Ecc/E"^""^ are 
considered [igj. but the theoretical results rely heavily on 
the phenomenological shape function. It is possible that 
the observed end point behavior of J/ip and the large 
ratio Rcc might indicate that the color-octet matrix el- 
ements are much smaller than previously expected. To 
test this thought we assume the color-octet contribution 
to be ignored and only consider the color-singlet contri- 
butions. Under this assumption, the ratio becomes 



Rc 



(j[J/^ + cc]/{a[J/^ -f cc] + a[J/^j + gg]). (2) 



In the following we concentrate on a{J /ip -\- gg) m 
NRQCD. Aside from the LO calculations in NRQCD (see 
related references in [3,111), Ref.;20] considered a[J/ip+ 
gg] within the framework of soft coUinear effective the- 
ory (SCET), and Ref. '^I'l summed over the leading and 
next-to-leading logarithms in the end point region of 
(t[J/'0 -f gg]. However, considering the crucially impor- 
tance of the NLO QCD corrections found in many heavy 
quarkonium production processes [1, 0, S, [HI [13, [IS], it 
is necessary to carry out the calculation of NLO QCD 
correction to e~'"e~ — > J/ij} -\- gg. 

We now present this calculation. We use FeynArts [2^ 
to generate Feynman diagrams and amplitudes, 
FeynCalc [i^] to handle amplitudes, and LoopTools 
to evaluate the infrared-finite scalar Passarino-Veltman 
integrals. Feynman diagrams for the Born, virtual 
correction, and real correction are shown in Fig. [1] 
Fig. [21 and Fig. [3l Note e+e~ J/^gcc is excluded 
in the real correction, because it should be included 
in the J/^ production associated with open charm 
e+e^ J/tfj + cc + X. Moreover, we include ghost 
diagrams in the real correction because we choose 
unphysical polarizations for the gluons in the final state. 

There are generally ultraviolet (UV), infrared (IR), 
and Coulomb singularities. Conventional Dimensional 
Regularization (CDR) with D = 4 — 2e is adopted to 
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FIG. 1: Three of the six Bom diagrams for e^e — » J/ipgg. 
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FIG. 2: Feynman diagrams for the virtual correction to 
e+e" ^ J/Tpgg- 



regularize them. 

The UV-divergences from self-energy and triangle di- 
agrams are removed by renormalization. The renormal- 
ization constants Zm, Z2, and Z3, which correspond re- 
spectively to the charm quark mass m, charm field ipc 
and gluon field A^, are defined in the on-mass-shell(OS) 
scheme, while Zg corresponding to the coupling as is de- 
fined in the modified-minimal-subtraction(A/5') scheme 
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where Nf 



(^)^r(i 



e) is a overall factor in our 



calculation, Poi^f) — — ^TpUf is the one-loop co- 

efficient of the QCD beta function, Uf = 4 is the number 
of active quark flavors, n;^ = 3 is the number of light 
quark flavors, and fi is the renormalization scale. 

IR singularities coming from loop-integration and 
phase space integration of real correction are found to 
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FIG. 3: Feynman diagrams for the real correction to e^e 



cancel each other. We use the method in [2g| to separate 
the soft and collinear singularities in the virtual correc- 
tions, and then treat the singular part analytically while 
the finite part numerically. 

We use phase space slicing method [l^l to extract poles 
in the real correction. The method introduces a soft cut 
Ss and a hard collinear cut Sc to the phase space. Then 
the cut region can be partly integrated, and becomes 
some color connected born cross sections multiplied by 
singular factors. While the remaining region, hard non- 
coUinear region, which is non-singular, can be integrated 
using the standard Monte-Carlo techniques. In order to 
make the method effective, Sc <C Ss is needed as men- 
tioned in 27|. With a careful treatment for the two cuts, 
we verified that our result is independent of the two cuts 
in a large range. 

We then find that by considering all NLO virtual and 
real corrections, and factoring the Coulomb singular term 
into the J ftp wave function, we get an UV and IR finite 
cross section for e+e" J/tp + gg. 

In the numerical calculation we use y/s — 10.6GeV, 
= 338MeV, and J/ijj wave function squared at the 

origin |i?j/^(0)p = l.OlGeV'^, which is extracted from 
the J/tp leptonic widthjil] at NLO in a^: |i?j/^(0)p = 



9Mi 



r6Q^[l-(16/3)Qs/7r] 

Taking Mj/^ = 2m (at LO in v) and to = 1.4 GcV, 
we get as(/i) = 0.267 for [i ~ 2rn, and the cross sec- 
tion at NLO in as is ^(e+e^ ^ J/ipgg) — 0.496 pb, 
which is a factor of 1.19 larger then the LO cross section 
0.418 pb. If we set ^ = ^/s/2, then Us =0.211, and the 
cross section is 0.394 pb. Since the experimental data 
correspond to the J/ij} prompt production, in addition 
to the direct production discussed above, we should also 
include the feeddown contributions from higher charmo- 
nium states which decay into J /ip- Since for the P-wave 
states Xcj the direct production rates in the non-cc as- 
sociated process are suppressed (e^'^e" ^ 7* ^ Xcjgg 
are forbidden due to charge parity conservation) , and for 
the il){nS){n > 2) states the decay branching ratios into 
J/V' + X are negligible, we only need to consider the 
'4>{2S) feeddown contribution. This implies an additional 



enhancing factor of 1.355 should be multiplied [ll|. In 
Fig. m we show the prompt production cross sections at 
LO and NLO as functions of the renormalization scale 
/i. We see that NLO QCD correction substantially re- 
duces the dependence, and enhances the cross section 
by about 20%. 

The NLO cross sections CT(e"'"e^ J/ipcc) were cal- 
culated in Ref. jll|. We list the values of prompt cross 
sections [m in TABLE HI together with the prompt cross 
sections (T(e"*"e~ — )■ J/iJjgg) obtained above. Then we can 
get the ratio Rcc at LO and NLO in . The dependence 
of Rcc on the renormalization scale ^ is shown in Fig. [5l 
where to = 1.4 GeV is fixed. The fi dependence for 



cr(e"'"e — > J/^gg) is mild, while for CT(e^ 



J/ipcc) 
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FIG. 4: Prompt cross sections of e"'"e~ J/ipgg as functions 
of the renormalization scale fi at LO and NLO in Os- The 
upper curves correspond to m = 1.4 GeV, and the lower ones 
correspond to m = 1.5 GeV. 





= 2.8GcV 


^l = 2.8GeV 


^l = 5.3GeV 


^l = 5.3GeV 




LO 


NLO 


LO 


NLO 




0.57 


0.67 


0.36 


0.53 


a{cc) 


0.38 


0.71 


0.24 


0.53 


Rcc 


0.40 


0.51 


0.40 


0.50 



TABLE I: Cross sections of prompt (feeddown included) 
J/ipgg (this Letter) and J/ipcc (Re.Tlj) production in e^e~ 
annihilation at B factories in units of pb. 



strong. A reasonable choice should be between /i — 2mc 
and jU = y/s/2, and more preferably the latter. 

Finally, we note that the large logarithms of log(l — z) 
appear at the endpoint in NLO calculation, where z — 
Ej/^/E'^/ ^^. The leading logarithms (LL) have been re- 



sumed in [20l |21| | . Using a similar approach, we define 



the differential cross section (and other quantities) as [21 1 
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rfo-rnsum I , (4) 

/ LO(NLO) 
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FIG. 6: The cross section of e^e~ — > J/tpgg as functions of the 
J/tp momentum Pj/^. Here n = 2.8 GeV and m = 1.4 GeV. 
The solid curve is the NLO+LL prediction, and the dotted, 
short-dashed, and long-dashed curves are the NLO, LO-I-LL, 
and LO predictions respectively. 



where (dcrrcsum)LO(AfLO) means expanding dcTrcsum in as 
to LO (NLO). To be consistent with our previous calcu- 
lation, we choose /ic — •\/2(l — z)i.ih, and jin = M ~ 2m. 
In Fig. [6] we show the cross sections of e^e~ J/'fpgg 
as functions of the J/ip momentum Pj/^. The correction 
of LL resummation to the total cross section is about 
—6.6% at LO and 0.5% at NLO. But it becomes large 
at the end point region when z — s- 1, suppressing the LO 
cross section and enhancing the NLO cross section. The 
LL resummation changes the J/i/i momentum distribu- 
tion near the end point, but has only a little effect on the 
total cross section. It is interesting to note that with the 
NLO correction, the J /ip momentum spectrum becomes 
much softer than the LO result. 

In summary, we find that by considering all NLO vir- 
tual and real corrections, and factoring the Coulomb sin- 
gular term into the cc bound state wave function, we get 
an ultraviolet, infrared, and coUinear finite cross section 
for the direct production of e^e~ —> J/ip + gg at y^s = 
10.6 GeV, which enhances the cross section by about 
20%. Adding the feeddown contribution from ijj{2S) the 



prompt production cross section of 
NLO in as is found to be (0.67:^°;^^) 
(0-53+o:?2) pb for ^i = with m 



"e J/tp + gg at 
pb for /i = 2m and 
_i_Q p — V / ^ w iLii III, = (1 .4 zb 0. 1) GeV. 

Together with the calculated CT(e+e~ J/^cc) at NLO 
we get Rcc « 0.50. The result significantly 



reduces the discrepancy between theory and experiment. 
The effect of the leading logarithm resummation near the 
end point on the J /ip -\- gg total cross section is found to 
be small. 

Note. Very recently Belle reported a new measurement 
with higher statistics [29|: 



a(e^ 



e- J/iP + cc) = (0.74 ± Qmtof^)ph, (5) 
J/tjj + non(cc)) = (0.43 ± 0.09 ± 0.09)p6.(6) 



FIG. 5: Rcc. as a function of the renormalization scale ^ at 
LO and NLO in as- Here rUc = 1.4 GeV. 



The observed cross section of e^e 



J/ip + non{cc) and 



Rcc are displayed in Fig. 0] and Fig. [5] with central values 



4 



and error bands in comparison with theoretical predic- 
tions. We see that, our predictions (NLO with feeddown) 
for a{e~^e~ — *■ J/^-\-gg) are consistent with the new mea- 
surement of a{e~^e~ — > J/ip + non(cc)) within certain 
uncertainties. Moreover, the predicted J/ip momentum 
spectrum also agrees with the experiment [2^ . Impor- 
tantly, our result of cr(e"'"e~ J/ip + non{cc)) indicates 
that the calculated cr(e+e^ J/tp+gg) has already satu- 
rated the observed cr(e+e^ — + J/'0-|-non(cc)), hence leav- 
ing little room for the color-octet contributions. These 
are also confirmed by a similar studvjsoj. which agrees 
with ours. 
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